Weighted Quadrature Rules for Finite Element Methods
نویسندگان
چکیده
We discuss the numerical integration of polynomials times exponential weighting functions arising from multiscale finite element computations. The new rules are more accurate than standard quadratures and are better suited to existing codes than formulas computed by symbolic integration. We test our approach in a multiscale finite element method for the 2D reaction-diffusion equation. Standard finite elements usually fail to accurately solve equations with multiscale behavior. This can happen if coefficients are oscillatory or if a small parameter multiplies some of the terms in the equation. A strategy to overcome these difficulties is to use special spaces instead of the space of piecewise polynomial functions [1, 2]. However, for polynomial basis functions, standard quadratures are exact and this is not the case for more complicated spaces. We investigate quadratures to integrate elementwise products of polynomials and exponential basis functions. Such integrals appear when developing enriched methods for reaction-advection-diffusion equations [2], but also in other contexts [3]. Quadrature formulas are simpler to implement into existing finite element codes than results of symbolic integrations. We define an N -point weighted quadrature in [a, b] with weighting function w by a set of integration weights Al and integration points xl ∈ [a, b] such that ∫ b
منابع مشابه
Applications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
In this paper Gauss-Radau and Gauss-Lobatto quadrature rules are presented to evaluate the rational integrals of the element matrix for a general quadrilateral. These integrals arise in finite element formulation for second order Partial Differential Equation via Galerkin weighted residual method in closed form. Convergence to the analytical solutions and efficiency are depicted by numerical ex...
متن کاملWeighted quadrature rules with binomial nodes
In this paper, a new class of a weighted quadrature rule is represented as -------------------------------------------- where is a weight function, are interpolation nodes, are the corresponding weight coefficients and denotes the error term. The general form of interpolation nodes are considered as that and we obtain the explicit expressions of the coefficients using the q-...
متن کاملGeneralized Gaussian Quadrature Rules in Enriched Finite Element Methods
In this paper, we present new Gaussian integration schemes for the efficient and accurate evaluation of weak form integrals that arise in enriched finite element methods. For discontinuous functions we present an algorithm for the construction of Gauss-like quadrature rules over arbitrarily-shaped elements without partitioning. In case of singular integrands, we introduce a new polar transforma...
متن کاملOptimally Blended Spectral-Finite Element Scheme for Wave Propagation and NonStandard Reduced Integration
We study the dispersion and dissipation of the numerical scheme obtained by taking a weighted averaging of the consistent (finite element) mass matrix and lumped (spectral element) mass matrix for the small wavenumber limit. We find and prove that for the optimum blending the resulting scheme (a) provides 2p+4 order accuracy for p-th order method (two orders more accurate compared with finite a...
متن کاملGeneralized Gaussian Quadrature Rules for Discontinuities and Crack Singularities in the Extended Finite Element Method
New Gaussian integration schemes are presented for the efficient and accurate evaluation of weak form integrals in the extended finite element method. For discontinuous functions, we construct Gauss-like quadrature rules over arbitrarily-shaped elements in two dimensions without the need for partitioning the finite element. A point elimination algorithm is used in the construction of the quadra...
متن کامل